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ON THE RELATIONSHIP BETWEEN SMALL-SCALE WIND SHEARS 
AND WIND PROFILE POWER SPECTRA 

SUMMARY 


Statistical wind shear information is used in the construction of 
synthetic wind profiles for the design of space vehicles. In practice, 
design scalar shears for various wind speed categories at various alti- 
tudes at which the shears are to be applied are derived from rawinsonde 
wind profile data. These data provide wind shear information for 
intervals equal to and greater than 1 km. To obtain wind shear infor- 
mation for intervals less than 1 km, detailed wind profile information 
like that obtained from the FPS-16 radar /Jimsphere wind sensing system 
is required. However, there are space vehicle launch sites which do 
not have detailed wind profile data in sufficient quantity to calculate 
statistical design shears. Accordingly, alternative procedures for 
deriving small-scale wind shears are required. The introduction of 
suitable hypo Lueses about the behavior of wind profile spectra at large 
wave numbers and the invariance of the distribution function of wind 
shears with regard to varying shear intervals facilitates the deriva- 
tion of an extrapolation procedure whereby one can extrapolate the 
rawinsonde shear statistics down into the region of shear intervals 
which are less than 1 km. To do this, it is assumed that any scalar 
wind profile from an ensemble of wind profiles can be represented in 
terms of a Fourier integral. The Fourier integrals permit us to 
relate the ensemble variance and mean of the wind shears to the power 
spectrum of the wind profile ensemble. The Fourier representation is 
general and accounts for the vertical nonnomogeneous statistical 
properties of the wind profile ensemble. It is hypothesized that (1) 
the power spectrum of the wind profile behaves like k" 2 * 4 at sufficiently 
large values of the vertical wave number k , (2) the vertical variation 
of the wind profile power spectrum can be neglected locally in the cal- 
culation of the ensemble mean and variance of the shear, and (3) the 
probability distribution function of the standardized shear variate is 
invariant with shear interval Az. These hypotheses are used to calcu- 
late the asymptotic behavior of wind shear statistics for sufficiently 
small values of Az. It is shown that the standard deviation and mean 
of the shear ensemble, as well as the shear for any percentile level 
occurrence, behave like (Az ) 0#7 for sufficiently small values of Az. 

This result is in excellent agreement with shear data obtained from 
FPS-16/Jimsphere wind profiles measured at Cape Kennedy, Florida. 


I. INTRODUCTION 


This report is concerned with establishing the relationship between 
small scale wind shears and detailed wind profile spectra in the first 
20 km of the atmosphere. Small scale wind shears are defined here to 
be vertical differences between the zonal, meridional, or scalar winds 
over vertical height intervals Az < 1 km. To this list of shears, one 
could add the vector shear, which is the magnitude of the wind shear 
vector whose components are the zonal and meridional wind shears. We 
will confine our discussion in this report lo scalar wind shears. 

Statistical wind shear information is used in the construction of 
synthetic wind profiles for the design of space vehicles. A detailed 
account of how a synthetic wind profile is constructed can be found in 
reference 1. Within NASA the practice has been to specify design wind 
shears based on scalar wind shear statistics? Moreover, two types of 
shears are considered: build-up and back-off shears, defined as 


build-up: 

SiCAZjZ) = V(z) - V(z 

- Z\z) 

(1) 

back-off: 

S 2 (Az,z) =V(z+Az) - 

V(z) , 

(2) 


where V(z) is the scalar wind at height z. The 99 percentile build-up 
and back-off shears enveloped over all altitudes are used for the design 
of space vehicles. To do this, the 99 percentile build-up and back-off 
shears are calculated from the annual empirical conditional scalar 
build-up and back-off shear distribution functions 


build-up: F x = F^SxCAz.z) |v x < V(z) < V 2 ) (3) 

back-off: F 2 = F 2 (S 2 (Az, z) | V x < V(z) < V 2 ) , (4) 

where and V 2 take on assigned values. Usually V 2 - = 10 m sec -1 

and 0 < Vi <100m sec -1 . This corresponds to nine wind speed categories 
at the reference level z. The distribution functions F x and F 2 are cal- 
culated from rawinsonde wind profile data and the smallest value of 
Az = 1 km. Thus, for example, the design scalar shears for Az ^ 1 km 
for the Eastern Test Range (ETR) are based on serially complete rawin- 
sonde wind profile data which consist of an eight-year sample of wind 
profiles spaced 12 hours apart with wind speed and direction specified 
at 1 km intervals in the vertical up to an altitude of 27 km observed 
during the period 1956 through 1963. The 99% values of S x (Az,z) and 


*These shears approximate the vector wind shears for sufficiently large 
wind speed magnitudes at the reference level z due to the relatively small 
directional variation in the wind vector with height at relatively high 
wind speeds . 
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s 2 (Az,z), S^AZjZ) and S 2 (Az,z) say, are calculated at 1 km height 
intervals between z = Az and z = 20 km for the build-up shears and at 
1 km height intervals between z = 0 and z = 20 km - Az for the back-off 
shears. The maximum value of Az is usually 10 km. This, for example, 
will result in 20 values of S+ (1 km, z) for Az = 1 km and 16 values of 
(5 km, z) for Az = 5 km. The design scalar shears are then taken to 
be the supremum values of S"![(Az,z) and S+(Az,z), sup (S^(Az , z) ) and 
sup (S 2 (Az, z)) , say. These calculations can be performed for various 
intervals (V x < V(z) < V 2 ) of the conditionalizing parameter V(z), 
resulting in two sets of design shear curves. One set will consist of 
the design build-up scalar shears, sup (S+(Az, z) ) , as functions of Az 
for various categories of V(z), and the other will consist of the design 
back-off scalar shears, sup (S^(Az, z) ) also as functions of Az for various 
categories of V(z). 

Because of the response properties of the rawinsonde system and the 
procedures used to calculate wind speeds from rawinsonde data, the cur- 
rently available rawinsonde data are not suitable for calculating design 
wind shear statistics for Az < 1 km. To calculate wind shear statistics 
for Az < 1 km. one requires detailed wind profile data like that obtained 
from the Jimsphere/FPS-16 radar wind sensing system. The number of 
detailed wind profiles from a particular launch site may not be large 
enough to calculate the required empirical conditional distribution 
functions, equations (3) and (4). However, it is possible to determine 
an extrapolation procedure from a smaller sample of detailed wind pro- 
files whereby one can extrapolate the values of sup (S*(Az, z)) and 
s up ( S 2 ( Az , z ) ) calculated with rawinsonde data down into the region 
Az < 1 km, and thus infer values of design scalar wind shear for vertical 
scales of distance Az < 1 km. 


II. ENSEMBLE STANDARD DEVIATION OF WIND SHEAR 


Let us consider an ensemble of scalar wind speed profiles. Each 
profile extends over the semi- inf inite domain 0 < z < oo. We select the 
ensemble of scalar wind profiles V(z) such that the scalar wind speed 
at z r , V(z r ) , in each profile lies in the interval V-l < V(z r ) < V 2 , 
where and V 2 have assigned values. We assume that we can express any 
wind profile from the ensemble with a Fourier integral 


V(z) 


00 

= J C(K)e' iKz d/c. 


-00 


( 5 ) 
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where the Fourier amplitude C (x) at wave number K (rad trT 1 ) is given by 


00 

C(K) v ( z ) elKZ dz - (6) 

-00 

The ensemble average of V(z) is 

00 

<V(z)> = J <C(K)>e" 1/<:z 6k, (7) 

-00 

where < > denotes the ensemble average operator. Upon subtracting (7) 
from (5), we obtain the scalar wind fluctuation v(z) at z with respect 
to the ensemble mean scalar wind profile; i.e # , 


v(z) 


00 

= f [C(K) - <C(K)>] e" iKz dK. 
-00 


( 8 ) 


Thus, any wind profile realization from the ensemble is given by 


V(z) = <V(z)> + v(z) . (9) 


By definition, the build-up and back-off scalar shears associated with 
the interval Az are 


build-up: 

S^ZjAz) =<S 1 (z,Az)> + 

s i(z, Az) 

(10) 

back-off: 

S 2 (z,Z^z) = <S 2 (z,Az)> + 

s 2 (z,Az) , 

(ID 


where 

<S 1 (z,Az)>= <V(z)> - <V(z-Az)> 
<S 2 (z,Az)> = <V(z+Az) - <V(z) > 
Si(z,te) = v(z) - v(z - Az) 
s 2 (z,Az) = v(z + Az) - v(z) 
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In the mathematical treatment that follows we will only consider the 
buildup scalar wind shears since we can formally obtain the back-off 
shear from the definition of the buildup shear by replacing z and z - Az 
with z + Az and z, respectively. Thus, any result we obtain for the 
buildup shear can be transformed into a statement about the back-off 
shear. 

We now proceed to calculate the variance of s x (z r ,Az). We can 
express the fluctuation of the scalar wind shear associated with the 
reference level z r in terms of the Fourier integral 

00 

s x (z r5 A 2 ) = J B(K?'')e~ ll< '' Zr [1 - e 1K,,Az ] d k", (13) 


where 


B(k) = C(/c) - <C(k)> . 


(14) 


We also can express in terms of the complex conjugate of the Fourier 
amplitudes C* (k) , so that 


SiCz^.Az) = J B *(/c) e XKZr [1 - e" 1KAz ] d K. (15) 


By definition <s !> = 0, so that the variance of s ± can be obtained by 
multiplying (13) and (15) and then performing an ensemble average. This 
yields 

00 00 

ol(z r ,£z) = J /<B*(k)B( 0> e^’> Zr (l-e 1K " Az )(l-e" 1 ' <Az )d/<"dK, 

-°0 -00 ( 16 ) 


where it is understood that V x < V(z r ) < V 2 . We can rewrite equation (16) 
by setting 


K - K" = K* , 


(17) 
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so that 


^(z Az) = r f<B*(K)B(K-K')>e lK ' Zr (l-e l(K ~ K ' )Az )(l-e' iKAZ ) d/c'd/c. 

1 C 


J 

00 -oo 


(18) 


Equation (18) gives the ensemble variance a 2_( z r >Az) °f the scalar wind 
shears as a function of the Fourier amplitudes of the scalar wind field. 

We will use (18) to obtain an approximate expression for a x (z r ,A 2 ) which 
will enable us to make estimates of the scalar shears for vertical scales 
Az < 1 km upon specifying the wind shears for Az H km for any percentile 
level of occurrence. This will become clear later. 

The power spectrum of v(z) at altitude z is 


0(K, z) = / <B^(/c)B(/c-K*) > e 


l/C z 


d K' 


(19) 


Although the quantity 0(K 9 z r ) is complex, the real and imaginary parts of 
0(k,z t ) 9 0 t (k 9 z r ), and 0±(k 9 z x ) are odd and even functions of K , respec- 
tively (see appendix A). Equation (19) permits us to express (18) in the 
form 


a^z^Az) = / 0(K i z)(l-e~ XKAZl ) d/c + / 0(k 9 z^Az) (1- e 1/<Az ) d k 9 


( 20 ) 


or in terms of the real and imaginary parts of 0(k 9 z x ) and 0(k 9 z r - Az): 


o-^z^Az) = 


[0 r (^,z r ) + 0 r (K, z r 


Az) ] [1 - cos (/cAz) ] 


+ [0 l (k, z^ - Az) - 0 i (K,z r )] sin(/cAz) f d/c 


+ i J' -^[0 ± (K, z r ) + 0 i (/c,z r - Az)][l - cos(kAz)] 

-00 

+ [0 r (K, z r ) - 0x(k, z r -Az)] sin(/<Az)T d/c. 


( 21 ) 



The imaginary part of this expression vanishes because the integrand of 
the imaginary part of (21) is an odd function of K on the interval 
-oo < z < oo # Thus, the reality of a^(z r ,Az) is guaranteed and 


a^(z r ,£z) 


00 

= / -f[0 r (K, z r ) + 0 r (K,z r 

-00 


- Az) ] [1 - COS (kAz) ] 


+ [0 i (K,z r - Az) 


0 i (*> z r )] 


s in (kAz) j>- d /c . 


( 22 ) 


If the ensemble of wind profiles has statistically homogeneous properties, 
then the spectrum of v(z) is real and independent of position z. This 
means that 


■KO = 0 r (K, z r ) = 0 r (K,Z r - Az) 
0 i ('<,z r ) = 0 i ( K > z r " Az) = 0 


(23) 


so that (22) reduces to 


a|(£z) = 4 / (1 - cos ( k £ z ) ) o(/c) d/c, 


(24) 


where we have used the even properties of $(k) to obtain an integral over 
the semi-infinite interval 0 < K < oo. The result given by (24) can also 
be obtained by noting that 0(k 9 z) can only be independent of z if 


<W' (k)R(k- K* ) > = ®(k) 8(k'). 


(25) 


Substitution of (25) into (18) yields the result given by (24). 
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Equation (22) expresses the general relationship between the ensemble 
variance of the wind shear and the spectrum of v(x). Equation (24) is a 
special case of (22) for the statistical homogeneous case and relates the 
ensemble variance of the wind shears to the homogeneous power spectrum. 

We will use these results to infer the dependence of a x ( z r ,Az) on Az 
from the power spectrum of the wind field. 


III. ENSEMBLE MEAN WIND SHEAR 


Let us now consider the ensemble mean shear. By definition, the 
ensemble buildup mean shear associated with the reference height z r and 
the shear interval Az is 


<Si(z r ,Az)> 


J <C(K")> e~ iK " Zr (1 


-00 


iK M Az 

e 


) d k". 


(26) 


Similarly, the ensemble mean shear can be expressed as 


<Si(z r ,Az) > = 


<C V) 


> e 


ixz r 


(1 - 


-i/CAZv , 

e ) d/c. 


(27) 


Multiplication of (26) and (27) will yield the square of the ensemble mean 
shear, namely, 


<S!(z r ,Az)> 2 = F (<C* M XC(,<-)>e^ r (1 - e i( -K- K ")te) 


-oo — oo 


• (1 


-iKAz 

e 


) d/c'dK. 


(28) 


Now the power spectrum of <V(z)> at altitude z is 


4r(K,z) 


OO 

= J' <C*(0><COk')> e iK ' Z d/c* . 
-00 


(29) 
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As in the case of the wind shear fluctuations, \|r(K,z) is complex, and the 
real and imaginary parts of this function, \[r r (K, z) and i|^(k, z), are even 
and odd functions of K, respectively. The proof of this statement is 
similar to the one in appendix A. Accordingly, the equation for 
<S 1 (z r ,Az)> 2 in terms of the real and imaginary parts of z) can be 

obtained from (22) by replacing 0 r *s and 0^ s in (22) with the corre- 
sponding \|r r f s and ty^'s, so that 


<S i(z r ,Az)> = 


|> r (K,z r ) +t r (*> z r - Az)][l - cos(kAz)] 


+ [y ± (K, z r - Az) - y ± (K, z r )] sin(/<Az) 


d/c. 


(30) 


In the case of statistically homogeneous conditions, the spectrum of 
<V(z)> is real and independent of z. Moreover, the spectrum is a delta 
function, so that 


ty r (^> z r ) — ~ Az) <V> &(/<) 

(31) 

^(k.z r> = = °» 

(32) 


where <V> 2 is the square of the ensemble mean wind profile which is 
independent of z. Substitution of (31) and (32) into (30) will show 
that <S 1 (z r , Az) > 2 vanishes for statistically homogeneous conditions. 

Equation (30) is the general relationship between the ensemble mean 
buildup scalar wind shear and the spectrum of the ensemble mean scalar 
wind profile. 


IV. STATISTICALLY HOMOGENEOUS WIND PROFILE ENSEMBLE 


We have found that, for an ensemble of statistically homogeneous 
wind profiles, the ensemble mean wind shear vanishes and the variance of 
the wind shear is given by (24), which we will repeat here for conven- 
ience: 


9 


oo 


(24) 


a^(Az) = 4 


[1 - cos (kAz)] $(k) d/c. 


o 


It is worthwhile to examine the implication of this formula for spectra 
that obey power laws. Power spectral calculations with detailed zonal, 
meridional, and scalar wind profiles observed with the FPS- 16/Jimsphere 
wind-sensing system show that, for sufficipntly large wave numbers, the 
power spectra of individual profiles behave like k" P, where p is a posi- 
tive number with numerical value approximately equal to 2.4 [1,2]. In 
these calculations the profiles were assumed to be statistically 
homogeneous. If we invoke the ergodic hypothesis (ensemble averages can 
be replaced with space averages), then, for sufficiently large wave 
numbers , 


4 >(*) = Pk" P , 


(33) 


where g is a positive quantity that depends on the intensity of the wind 
fluctuations about the ensemble mean wind profile, which in this case is 
independent of z . 

We wish to calculate the ensemble variance of the shears for small 
values of Az, Az < 1 km, say. Now for sufficiently small values of Az , 
the behavior of Q(/c) at small wave numbers is relatively unimportant in 
the evaluation of (24) because the factor 1 - cos ( kAz) in the integrand 
tends to suppress the contributions to a^(Az) from ®(k). Accordingly, 
we will use (33) to evaluate (24) over the complete interval 0 < K < oo, 
so that 


o^(Z\z) 


4p 


K _P [ 1 


- cos(kAz)] d/c. 


or 


Oi(Az) 


(3 2jt cosec (y jt(p-l)) 

(p-i) r(p-i) 


(Az) 


P-1 


(34) 


(35) 


provided 1 < p <3, which corresponds to the situation we are dealing with 
here. Equation (35) shows that, for statistically homogeneous conditions, 
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the ensemble standard deviation of the scalar wind shear depends on £z, 
asymptotically, like 

Ell 

<*i “ 2 , (36) 


where 


'1,0 


n cosec ittp-l))-^ 1 / 2 


(P-1) r(p-i) 


^/ 2 . 


(37) 


V. NONHOMOGENEOUS WIND PROFILE ENSEMBLE 


In reality wind profiles are not homogeneous* One example of the 
nonhomogeneous behavior of wind profiles is the general increase and 
decrease of the scalar wind speed below and above the core of an atmos- 
pheric jet stream* A second example is the change. in the eddy structure 
as height increases through the tropopause; i.e., the wind fluctuations 
about the running mean of a wind profile realization above the tropopause 
tend to be larger than the ones below the tropopause, especially when the 
troposphere is convectively unstable. Gravity wave theory provides a 
third example. According to Hines [3], the Fourier components 


u(k,x,t) 


of the zonal and meridional wind at position x, and time t in a gravity 
wave are given by 



U(k, W (k))e z/2H 


(38) 


and 
is 

the Fourier amplitude. Hines [3] obtained this solution of the inviscid 
hydrodynamic equations by perturbing a quiescent isothermal atmosphere 
with infinitesimal perturbations. The important point here is that Hines’ 
gravity waves constitute one plausible class of motions in the atmosphere 
which have Fourier amplitudes that increase with altitude like e 2 ^**, thus 
rendering the process nonhomogeneous. There are other types of plausible 
atmospheric wave motions that have Fourier components that vary with 


where k ik a wave number vector, H is the atmospheric scale height 


oj(k) is the eigen frequency associated with k. The quantity Ue z 


m 
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height in a manner similar to Hines' gravity waves. In view of the 
obvious nonhomogen eous nature of atmospheric wind profiles, it is reason- 
able to inquire about how the variance and mean of the scalar wind shear 
depends on Az. Unfortunately, the vertical variation of the power spectra 
of wind profiles is not well known. However, by making assumptions about 
how 0(K,z) and \|t(k,z) depend on K and z, we can make inferences about how 
a;L(z r ,Az) and <S 1 (z r ,Az) >, which we shall denote with depend on Az. 
Furthermore, the introduction of an assumption concerning the behavior of 
the probability density function of S-^ZjAz) as Az decreases will permit 
us to predict how S^z.Az) for any percentile depends on Az. The result- 
ing predictions can be tested with existing experimental data on the dis- 
tribution function of wind shears from Cape Kennedy, Florida. 

The integrands of (22) and (30) are even functions of K, so that we 
can replace the integrals over the interval -<» < k < °o with integrals 
over the interval 0 < K < thus , 


00 

of(z r ,Az) = 2 /{ [V K,z r 


) + 0 r (K,z r 


- Az)][l - cos(kAz)] 


+ [0 i (<,z r - Az) - 0 ± (K,z r )] sin(KAz)j- d k. 


(39) 


oo 

i^(z r ,Az) = 2 /{[ tr (,c.z r ) + \|r r (K,z r - Az) ] [1 - cos(kAz)] 


+ [i|f i (K,z r - Az) - \lf i (/C,z r )] sin(KZ^)j- d/<. 


(40) 


As stated previously, we do not know how 0(k 9 z) and \[r(/c,z) depend on z. 
However, we do know that these quantities are functions of altitude, and 
this means that there will be different spectra or sorting by wave number 
for each value of z. Nevertheless, the process described is a physical 
one, and we hope that the variation of the spectra from position to posi- 
tion will not be too abrupt, but rather even with smooth variation. With 
this philosophy in mind, we expand 0 r (/<,z r - Az), 0i(K,z r - Az), 

^r(K> z r “ Az) and \|ri(/C, z r - /\z) in Taylor’s series about the point z = z r 
for small values of Az so that second- and higher-order terms can be 
neglected. This permits us to express (39) and (40) as 
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oo 


J^(z r ,Az) = 4 / 0 r (K,z r ) 


1 - 


4 ( K > z r ) 

2 ^ r (/<, Zr ) 


Az 


[1 - cos(kAz)] 


0 [(K, z r ) 


20 


r (^ 7 T) ^ sin ^); dK > 


( 41 ) 


|i i( z r , ^ z) = 4 / *r* K,z r* 


1 - 


t^» z r ) 

2 t r (K> z r ) 


Az 


[1 - cos (/<Az) ] 


Z A 

2\|r r (/c,z r ) 


Az sin(/cAz) r d/c, 


(42) 


where prime denotes differentiation with respect to z. 
If we now assume that 


0 1 (K, z ) 

r ’ r 


20 r (K,z r ) 


Z^Z 


« 1 


(43) 


^(^,z r ) 


20 r (K,z r ) 


Az 


« 1 


(44) 


*r^ ,z A 

2 i r (<,z j .) 


Z^z 


« 1 


(45) 


i-(^, z r ) 

2 t r ('<» z r ) 


Z^z 


« 1 . 


(46) 
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We can then approximate (41) and (42) as 


oo 

o^(z r ,Az) = 4 J' 0 r (/c, z r )(l - cos(kAz)) d/c, (47) 

o 

oo 

>Az) = 4 y' t r ('<,z r )(l - cos(kAz)) dK. (48) 

o 

The result given by (47) is similar in form to the one for the statis- 
tically homogeneous case; however, (47) permits 0 r (/c,z r ) to vary with z r , 
but the variation must be sufficiently slow such that (43) and (44) are 
satisfied. Note that, in the nonhomogeneous case, p.-L(z r ,Az) is a func- 
tion of z r and Az, while in the homogeneous case, |i 1 is equal to zero. 

We now assume that, for sufficiently large wave numbers, 0r(K>z r ) 
and ^ r (K,z r ) depend on k through power laws of the form 


-p (z) 

0 X (K, z) = P r (z )k r 


(49) 


-q r (z) 

t r (K,z) = 7 r (z ) K , 


(50) 


where fa, 7 r , p r , and q r are positive definite functions of z. Substitu- 
tion of (49) and (50) into (47) and (48) yields 


3 r (z r ) 2 it cosec (^ «(p r (z r )-l)) p r (z r )-l 

* (MO-i) r( P )-D (to) 


(51) 


r r 


r r 


2 7 r (z r ) 2n cosec ( 2 (z )“!)) q r (z r )-l 

l a i^ Z r , ^ z ^ - ~ 7Z t.n "nTT, ~ (Az) 


(q r (z r )-i) r(q r (z r )-i) 


These results show that a x <x (Az) ^ Pr_1 ^ 2 and 
f iciently small values of Az. 


(52) 

Hx K (Az) ^^ 2 for suf- 
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Various investigators have experimentally examined the nature of the 
dependencies of and on Az. Adelfang, Ashburn, and Court [4] have 
shown that the annual climatological mean and standard deviations of the 
vector shear at Cape Kennedy, Florida, depend on the reference level z r 
and Az. In particular, they found that both the mean and standard devia- 
tions are proportional to (Az) 2 / 3 and the coefficients of proportionality 
are functions of height. Their analysis was based on 1194 Cape Kennedy 
Jimsphere profiles. The computations were made at z r = 8, 12 and 16 km 
for Az = 50 to 5000 m. Other authors have deduced power law relation- 
ships between Az and and a x (see, for example, references 5 and 6). 

In view of the results of Adelfang, Ashburn, and Count, we will set 


P r (z r ) = q r (z r ) = a. 


(53) 


where a is a constant. As mentioned previously in Section IV, power 
spectra of individual scalar wind profiles behave like at sufficiently 
large wave numbers, where p has numerical value approximately equal to 2.4. 
If we assume that this value of p obtained from individual profiles can 
be identified with a, thena - 1 = 1.4. Thus, we hypothesize that 


1 i 1 (z r ,Az) a ( Az ) ° • 7 
a^z-^Az) « (Az)°* 7 , 


(54) 


where the coefficients of proportionality are functions of z^. 

Let us now assume that the conditional probability density function 
of Sjl can be transformed into the conditional probability density func- 
tion: 


p = p(x | Vi < V(z r ) < V 2 ) , 


(55) 


where x is the nondimens ional variable: 


S 1 (z r ,A?) - m(z r ,Az) 

0 1 (z r ,Az) 


(56) 
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Integration of (56) over the interval -<» < x < x + will yield probability P 
that x § x + , given that < V(z r ) g V 2 ; i.e., 


P(x g x + | V x < V(z r ) < V 2 ) 



Vi < V(z r ) g V 2 ) dx. 


— 00 


Now 


(57) 


P(x g x + | Vi < V(z r ) S V 2 ) = P (S ! g St | V x < V(z r ) g V 2 ), 


(58) 


where S* is the value of S x associated with x + , so that the value of S* 
for an assigned value of P can be obtained by solving (56) for s\ with 
x = x + ; i. e. , 


st(z r> Az) = m(z r ,Az) + x + a 1 (z r ,Az). (59) 

If we now assume that the distribution function of x is invariant with 
Az for a fixed value of z r , then it follows that x + is only a function of 
P, and we conclude from (54) and (59) that 

sl(z r ,Az) = G(z r ,P(x g x + | Vi < V(z r ) S V 2 ))(Az)°‘ 7 , (60) 


where G is a function of z r and P. To obtain the design buildup wind 
shear, we select the supremum value of si[(z r> Az) for P = 0 # 99 from the 
set of 99 percentile buildup shears that is obtained by letting z r take 
on all permissible values. In view of the hypothesis that G is independent 
of Az, it follows that sup (S^(z r , Az)) is proportional to (Az)°* 7 . The 
results given by (51), (52), and (60) also hold for the back-off shears. 
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VI. EXPERIMENTAL DATA 


To verify the theoretical speculations and hypotheses presented in 
the previous sections, Mr. C. Hill of the Aerospace Environment Division 
has made available his calculations of the distribution functions of 
scalar wind buildup and back-off shears for Cape Kennedy, Florida. Hill 
performed his calculations with a sample of approximately 256 FPS-16 
radar/Jimsphere profiles observed during the period January 3, 1967 to 

June 30, 1967. The reference levels for his calculations were z r = 1,2, 

... 15 km. At each of these reference levels, he stratified the pro- 
files into 10 m sec " 1 wind speed categories such that V 2 - V x = 10 m sec " 1 
and V 2 - 10, 20, 30, .... He then calculated the buildup and back-off 
shears for each wind speed category for each value of z r = 1 , 2 , ... 15 km. 

Generally, his results show that the buildup and back-off shears for 
any particular percentile are functions of altitude, and increase like 
(Az)°* 7 for sufficiently small values of Az; i.e., Az < 1 - 2 km. The 

results for z r = 10 km for 10 < V(z r ) ^ 20 m sec " 1 and 30 < V(z r ) g 40 m 

sec " 1 were selected for presentation here, and they are typical examples. 
Figures 1 and 2 contain the buildup shears, and figures 3 and 4 contain 
the corresponding back-off shears. The figures show that the shear for 
a given percentile varies like (Az)° -7 for Az ^ 1 km. These results seem 
to verify the speculations and hypotheses presented in the previous sec- 
tions: (1) The ensemble fluctuation spectrum and the spectrum of the 

ensemble mean wind profile are sufficiently slowly varying functions, 
such that ]2 r (K,z r - Az) ~ 0 r (K,z r ), etc. for sufficiently small values 
of Az, (2) the real parts of the power spectra depend on wave number like 
K“ 2 - 4 , and (3) the distribution function of the standardized shears 
(Si - p 1 )/a 1 and (S 2 - P 2)/ a 2 are invariant with Az. 

The fact that S+ and behave like (Az)° -7 for sufficiently small 
values of Az permits the estimation of small scale wind shear for Az < 1 km 
with rawinsonde wind shear data at Az = 1 km. The procedure consists of 
specifying Az and calculating the associated shear with the formula 



9 Az) 


S^ > 2 (z r ,l km)(Z^)°' 7 , 


( 61 ) 


where £z is in units of kilometers. 
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100 1000 

Az, m 

Figure 1. Cape Kennedy, Florida buildup scalar shears S^(z r ,Az) as func 

tions of Az for z r = 10 km and 10 < V(z r ) i 20 m sec -1 for 
various percentiles . 
















VII . CONCLUDING COMMENTS 


The purpose of this report has been to determine the relationship 
between scalar wind shears and power spectra of atmospheric flows. We 
have seen that there is an intimate relationship between the ensemble 
mean shear and the spectrum of the ensemble mean wind profile and 
similarly between the ensemble standard deviation of wind shear and the 
spectrum of the wind fluctuations relative to the ensemble mean wind 
profile. The analysis accounted for the nonhomogeneous aspects of wind 
shear statistics, and in this respect the analysis is general. The 
hypotheses listed in Section VI enabled us to make relatively accurate 
predictions about how the buildup and back-off shears vary with Az for 
a given percentile level of occurrence. In particular, we found that 
the scalar wind shears depend on Az like (Az )°* 7 for sufficiently small 
values of Az. The analysis in this report can be applied to zonal and 
meridional wind shears by making slight changes in notation. 
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APPENDIX A 


Proof that the 0 v (k,z) and 0^(k, z) are Even and Odd Functions of k 


We assume that v(z) can be expressed as 


v(z) = / B (K 2 )e 


“1/Cp2 


dz. 


(A-l) 


where 


B(k 2 ) / v(z) e 1/<2Z dz 


(A- 2) 


and k 2 is a wave number. The complex conjugate of B (k) evaluated at 
wave number k x is 




J "* v ( ^ ) e d|. 


“00 


(A- 3) 


Multiplication of (A-2) and (A-3) and ensemble averaging the resulting 
relationship yields 


<B'' ( X i) B ( x 2 ) > - 


( 2 re) " 


<v(z)v(£)> e^^ 22 dz d| 4 


-00 -00 


(A- 4) 


Complex conjugation of (A-4) yields 


<B /V (k 1 )B(/<2) > — 


(27t)‘ 


< v(z)v(^)>e- i( ' < ^-'<i^). 


(A- 5) 
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The right side of (A-5) is <B*(-/C 1 )B(-K 2 )>. Therefore, we have 

<B*(k 1 )B(k 2 ) > = <B* (-/<!> B(-#C 2 ) > . (A- 6) 

Complex conjugation of (A-6) yields 


<B*(k x )B(/< 2 )> = <B*(-K 1 )B(/< 2 )>*. 


(A- 7) 


Now the power spectrum of v(z) at wave number k and altitude z is 


0(/c, z) = / <B*(k)B(/c - K’)>e iK ' Z d k* . 


(A- 8) 


The corresponding spectrum at wave number -/< is 


0(-K, z) = f<B*(-K)B(-K - «')>e lK ' Z d#c* 


(A- 9) 


Complex conjugation of (A-9) and utilization of (A-7) yields 


0 *(-< 


00 

,z) = jf<B* 


(/c)B (kH-zc* ) > e 


i \ ^ -iK 1 2 j i 

' “ d/c' . 


(A- 10) 


Let 


K 1 = -K". 


(A— 11) 
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Thus, (A- 10) can be written as 


0*(-K, Z) 


00 

= J<B*(k)B(k-k")> e 1K " z dK". 

-00 


(A- 12) 


However, k" is a dummy variable of integration which we can replace with 
K* in (A-12), so that 


/(-«, z) 


00 

= J<■E>*(K)B(K-K')>e iK,Z <lK•. 

-00 


(A- 13) 


The right-hand sides of equations (A-8) and (A-13) are equal. This means 
that 


0(k,z) = 0*(-K,z), 


from which we conclude 


0 r (K,z) = 0 t (- k , z) 

0 i (K,z) =-0 i (~K, z). 


(A- 14) 


(A- 15) 
(A- 16) 


Equations (A-15) and (A-16) state that the real and imaginary parts of 
0, 0j-, and 0j_ are even and odd functions of K, respectively. 
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